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Abstract 

We consider a two parameter family of Z2 gauge theories on a lattice discretiza- 
' tion T(A4) of a 3-manifold M. and its relation to topological field theories. Familiar 

£H ■ models such as the spin-gauge model are curves on a parameter space V. We show 

that there is a region To C T where the partition function and the expectation value 
(W j? (7)} of the Wilson loop can be exactly computed. Depending on the point of 
To, the model behaves as topological or quasi-topological. The partition function 
is, up to a scaling factor, a topological number of AA. The Wilson loop on the 
other hand, does not depend on the topology of 7. However, for a subset of To, 
(Wail)) depends on the size of 7 and follows a discrete version of an area law. At 
the zero temperature limit, the spin- gauge model approaches the topological and the 
quasi-topological regions depending on the sign of the coupling constant. 

^ 1 Introduction 

00 

I/"") ' A lattice gauge theory with gauge group Z2 is the simplest example of a gauge theory 

[I]. In dimension d = 2, the partition function for Z2 (as for any other compact gauge 
group) can be computed in various ways. In dimensions larger than two, however, 
the simplicity of Z2 does not help us to solve the model. Even without matter, the 
relevant models are non trivial and exact solutions are not known. Such solutions 
would be a very important achievement. A Z2 gauge theory on a cubic lattice can be 
made dual to the 3D Ising model, an outstanding problem in statistical mechanics 
[21 [3]. It is clear that in d = 4 the problem is at least as difficult as in d = 3. In this 
paper we will not have much to say about d = 4 since the tools we use are peculiar 
to dimension three. 

However difficult, 3D lattice models with local gauge symmetry are not always 
beyond the reach of exact solutions. That depends on the dynamics, i.e., the choice 
of an action for a lattice plaquette. Topological Quantum Field Theories (TQFTs) 
are examples where one can perform exact computations. Examples have been 
constructed on the lattice in dimensions d = 3 [U [5] as well as d = 4 [6] . Once the 
topology on the manifold in question is fixed, the partition function does not depend 
on the lattice size and can be trivially computed for a discretization with very small 
number of sites, links and plaquettes. Such models are very simple from the physical 
point of view. It follows from topological invariance that transfer matrices are trivial. 
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Despite being trivial dynamically for a fixed topology, TQFTs are quite relevant 
in physics. The reason being that TQFTs can come out as limits of ordinary field 
theories in the continuum as well as in the lattice. The most celebrated example is 
topological order in condensed matter physics [7] where the physics at large scales is 
described by a TQFT. Something of similar nature also happens for lattice theories 
in d = 2 that are quasi-topological [8l |9]. They are reduced to a TQFT at the 
appropriate limits. 

In order to understand the relationship between fully dynamical d = 2 models and 
their possible topological limits one can first look at quasi-topological models. They 
are very easy to work with since we can compute all relevant quantities. Quasi- 
topological models in d = 2 are a nice set of toy models for this purpose. In 
particular, the relation between the original models and their topological limits is 
made very explicit. As for d = 3 the situation is more complicated. We do no 
have at our disposal toy models that are at the same time not topological and easily 
computable. We have no choice but to work with fully dynamical theories where no 
exact computations are available. 

The focus of this paper is to investigate how lattice field theories with local 
7j2 gauge symmetry are related to topological theories. Before going any further 
we need to say what we mean by a lattice model being topological. Let T(M) 
be a lattice triangulation of a fixed compact 3-manifold A4. We say that such 
a model is topological if the partition function is the same for all triangulations 
T(A4). Actually, this is a weak definition since we may ask that not only the 
partition function but the expectation value of all observables to be of a topological 
nature. In any case, we have to go beyond the usual regular cubic lattices and take 
into account arbitrary lattices. 

In |10| we investigated the spin-gauge model in d = 3. We showed that in the 
limit /3 — > oo the partition function Z is given by 

Z{T{M)) |^oo = Z to P(M) 2 n *+ n l- n t, (1) 

where Nt,Nf and Nl are the number of tetrahedra, faces and links of T(M) and 
Z top (A4) is a topological number and as such does not depend on the discretization 
T(A4). Equation ([1]) tell us that at the limit f3 — > oo the partition function is not 
strictly speaking topological since it depends on the triangulation. However it does 
not depend on the details of the triangulation but only on its size. For this reason we 
say that the the partition function is quasi-topological. It follows from ([TJ that the 
partition function can be computed for all triangulations. Let Tq(M) be a lattice 
triangulation where the numbers Nt , Np and Nl are very small such that Z(Tq) 
can be written down explicitly. For an arbitrary lattice T(A4) we have 

Z{T)\^ = Zpo)^ 2 (N F -N Fo)+ (N L -N Lo) -(N T -N To) _ (2) 

We will find it convenient to rewrite Z as a product of local Boltzmann weights, 
in the form 

Z=z2Hw(f), (3) 

{9a} f 

where the produt is over all faces of the triangulation and sum is over all configu- 
rations. Let a, b and c be the links of a face / and (g a ,9b>9c) a g au g e configuration 
at /. The corresponding local Boltzmann weight for the spin-gauge model can be 
written as 

W^(g a ,g b ,g c ) = e^^. (4) 
Another very common choice is to set the local Boltzmann weight to 

W^(g a ,g b ,g c ) = e-^-^\ (5) 
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which corresponds to the usual gauge theories where flat holonomies will have the 
highest weight. 

As we will see in this paper, the relationship between gauge models and TQFTs 
can be better understood if we depart from a specific example as in [10J and consider 
a more general class of gauge models. In order to have a gauge theory, the local 
weight should depend only on the product of the gauge variables around an oriented 
plaquette: 

W(g a ,gb, 9c) = M(g a g b g c ). (6) 

Gauge invariance means that M(g) is a class function or, in other words M{hgh~ l ) = 
M(g), V7i G G. The character expansion for the group Z2 is very simple and implies 
that 

M(g a , 9b, 9c) = mi (g a gbg c ) + rn , (7) 

where mi 6 M. The original spin-gauge model with one parameter f3 can be recovered 
by restricting the model to a curve (mo, mi) = (cosh , sinh (/3)) in this two- 
dimensional parameter space. We will refer to the parameter space as V. 

The first question to be addressed is the generalization of equation ((2|). We will 
show that there is a subset Tq of the parameter space T such that the partition 
function Z(m,Q, m\, T(A4)) can be written as a product of a topological invariant 
Z top (A4) times a known function depending on the numbers Nt, Np and Nl of 
tetrahedra, faces and links. Again, an easy consequence is that at Tq the partition 
function Z(rriQ, mi, T{M)) can be computed for any lattice T(A4). The subset Tq 
is made of two pairs of lines, namely, rrii = 0, % = 0, 1 and mi = ±mo- It turns out 
that these two regions of Tq have different properties. For instance, the topological 
invariant that appears in the first pair is trivial and Z top (A4) = 1 for all compact 
manifolds Ai. As for second pair, Z top (A4) depends on the first group of co- homology 
of Ai [10J. The two pairs of solutions are also related to high and low temperature 
limits as it will be clear from the discussion on section [2j 

As in any gauge theory, one would be interested in more observables than just the 
partition function. In particular it is important to calculate the expectation value 
(W if, (7)) °f Wilson loops for arbitrary representations R of the gauge group and 
closed curves 7. In our previous work [10] we considered only the partition function. 
In the present paper, we would like to go further and ask whether (^(7)) can 
be computed for some points of V. As it happens for the partition function, such 
computation can be performed for all points of Tq. Note that in a truly topological 
gauge theory such as Chern-Simons, {Wr(^)) is a topological invariant of 7. That 
is not true for all points of Tq. It turns out that (Wr^)) depends on the size of 7 
for points of Tq of the form (mo, mi) = (A, —A). Since there is no dependence on 
the details of 7 we say that the observables (Wr(^)) are quasi-topological. 

Our approach is based on the fact that a large class of lattice models, topological 
or otherwise, can be described by a set of algebraic data on a vector space V. These 
data comprises of a multiplication m, a co- multiplication A and an endomorphism 
S such that S 2 = 1. It is also assumed that there is a unity e and a co-unity e. It has 
been shown in [H[5] that when the data (m, A, S, e, e) defines a Hopf algebra, one can 
construct a lattice topological field theory. We observed in |10| that the same data 
can be used to describe an ordinary Z2 gauge theory. In |10| . however, (m, A, S, e, e) 
is not a Hopf algebra. In particular, the co-multiplication is not an algebra morphism 
as it happens for Hopf algebras. Another important difference is that instead of a 
fixed algebraic data, we had a one parameter family m(/3) of multiplications where /3 
is the coupling constant of the model. It turns out that a Hopf algebra is recovered 
in the limit (5 — > 00 and the model becomes quasi-topological. It is also possible to 
have examples with matter fields via a family of co- multiplications A (A) where A is 



3 



the corresponding coupling constant In this paper, however, we will be limited 
to pure gauge theories. As stated before, we will consider gauge theories with two 
coupling constants given by ([7]). If we were to follow the formalism of [10J, that 
would be encoded in a two parameter family of multiplications m{mQ ) m\). The 
model considered in [5] corresponds to the unique point (mo, mi) = (1,1) where 
m(l, 1) together with A,S,e, and e define a Hopf algebra. For the present paper, 
however, the Hopf structure is less important. What matter are the points where 
the model is quasi-topological. That happens for (mo, mi) belonging to the region 
To described above. The data m(mQ, mi), A, S, e, defines a Hopf algebra only at 
a single point of Tq. Just as the topological model of [U [5], the limit f3 — > oo of 
equation ([JJ also corresponds to a point in the set Tq. 

The organization of the paper goes as follows: On Section [2] we explain how the 
algebraic data can be used to encode the model and how familiar models, such as 
the spin-gauge model, fit into the parameter space T. On Section[3]we determine the 
subset To where the model is quasi-topological. The computation of the expectation 
value for Wilson loops is investigated on Section UJ We show that (Wr(j)) can 
be computed for IV The model is not topological for all Tq. We show that for 
a particular region of To, {Wr(j)} depends on the size of 7 and follows a discrete 
version of an area law. 

We close the paper with some final remarks on Section [5l 

2 The Partition Function 

In this section we introduce the formalism we will use to describe the partition 
function of a Z2 gauge theories. We will loosely follow [U[5j[IU] making the necessary 
modifications to suit our purpose. 

Let T(A4) be a triangulation of a compact 3-dimensional manifold A4. For the 
description of a pure gauge theory, what is relevant in a discretization is the set 
of faces and how they are interconnected. To encode the connectivity information 
we will split T{M) into individual faces and record the information on how they 
should be put back together. This process can be described as follows. For each face 
fk £ T{M) we associate a disjoint face and for each link lj € T{M) we associate 
a hinge object Hj with rij flaps as illustrated in figure [TJ The number nj of flaps is 
equal to the number of faces of T(Ai) that share the link h. To reconstruct 

T{M) from the disjoint faces Fj., we can use the hinges Hj to determine which 
faces are to be joint together. This is illustrated by figure [TJ The faces F\. and the 
hinges Hj are to be given an extra structure called orientation. For the case of Z2 
gauge theories, this orientations are not relevant but we will mention them as they 
help to organize the model. We will call the set {F^,Hj} and its interconnections a 
decomposition of T(A4). 

We will use the decomposition {Fk,H m }, plus some extra data, to define a par- 
tition function. The first step is to choose a vector space V of dimension n. The 
edges of a face F^ have to be enumerated from 1 to 3. That amounts to a choice 
of orientation of the face and a choice of the starting point (see figure [2] (a)). The 
edges of l m G F k carry configurations ( 3), with £ {1, • • • , n}. A statistical 

weight M„k„k n k will be associated to the face Fu. Note that M„k h k^k can be viewed 
as the components of a tensor M S V <S> V <g> V. Furthermore, it should be invari- 
ant by cyclic permutations of its indices since we do not care which edge is to be 
numbered as the first one. On the other hand, a change in orientation can affect the 
corresponding weight since M a b c may not be the same as M c b a . In a similar fashion, 
the flaps of a hinge H m can be cyclically numbered from 1 to q = I(l m ). Once more, 
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(a) 



(b) 



Figure 1: In (a) we have a small piece of a triangulation. Figure (b) shows its corre- 
sponding decomposition into flaps and faces. 

this is equivalent to give H m an orientation, as illustrated by figure [2] (b). The flaps 
of H m carry configurations (a™, • • • , a™)> a™ ^ {■"■> '*' > n i J us ^ n ^e ^ ne edges of a 
face. That will correspond to a statistical weight A"! 1 '""™. 

Such numbers can be interpreted as the components of a tensor A € V* (g> • • • ® V * 
and have to be invariant by cyclic permutations of the indices. As before, change in 
orientation of H m will change the statistical weight to A a « %-i"' a i . 



Figure 2: The links on a triangular face are numbered from 1 to 3. The orientation is 
indicated in (a). Figure (b) shows a hinge with four flaps and its orientation. 
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(a) 



(b) 



Once we fix an orientation for each Fk and an orientation for each H m , we 
produce a tensor M aia2a3 (Fk) for each face and a tensor A bl "' bq (H m ),q = I(l m ) for 
each hinge. One can see that the product 



U^alWU^'^i 11 ^ (8) 



has one covariant index for each edge and one contra-variant index for each flap 
of {Fk,H m }. The partition function will be the scalar constructed by contracting 
all indices. A covariant index ah is to be contracted to a contra- variant index fj 1 



5 



whenever the corresponding edge Fk and flap H m are to be glued together. In other 
words, we define the scalar Z as 

z=U M 444mU AbT ''' bT ^ II (9) 

k m bt={gluings} 

where the last product is responsible for contracting indices. There will be a 5J (b) 

for each paring edge-flap (F r ,H s ) that are glued together, as illustrated by figure[3j 
We can simplify the notation of Q by eliminating the Kronecker deltas and writing 

z= n M abc (f) n a 6i - 6 ^(o, (io) 

/e{F} ie{L} 

where {F} and {L} are the set of faces and links of the triangulation. A contraction 
on the indices corresponding to gluings is understood. 




Figure 3: In terms of tensor, we interpret the gluings as contractions of the corresponding 
indices. 



As of now, the partition function ([9]) is not very useful. We need to be more 
precise about the weights M aia2CL3 (Fk) and A bl "' bn (H m ) if we want Z to be related 
to physical models like the spin-gauge model. Note that Z depends on the choice 
of orientations of individual faces and hinges. This dependence on the orientation 
should not be present in the final model. Furthermore, the weight function should 
be the same for all faces F^ and hinges H m . To go any further we need to constraint 
the tensors M a b c (Ffc) and A bl '" bn (H m ). That can be done with the help of some 
algebraic data that we will now introduce. 

The first algebraic structure we need is a product on V defined by 

cj) a ■ <j) h = M ab c <j) c . (11) 

where we are using a tensorial notation with the usual convention of sum over re- 
peated indices. We will use the symbol A to refer to the vector space V to emphasize 
that we are now working with an algebra. 

In this paper, we will choose A to be the group algebra of Z2. The group elements 
are written as ( — l) a , a = 0, 1 and a basis for A is {</>o, <j)\}. The product is defined 
by 

4>a- 4>b = <Pa+b- (12) 

Whenever a sum a + b of indices appear, as in (fT2|) . it will always denote sum module 
2. This product can also be given in terms of the tensor M ab c as 

M ab c = S(a + b,c). (13) 

It is also convenient to define the dual vector space A*, the dual base {if} 1 } and 
the usual pairing 

(rAb) = 5 a h . (14) 
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We then define the trace T € A* as 

T = M ab b r- (15) 
Given a face F^, as in figure [3j we define the associated weight M aia , 2a:j as 

M ai a 2 a 3 ( z ) = (T^a! ' 4>a 2 " </>a 3 " «), (16) 

where z = a°0o + a 1 '/'! is a generic element of A. Since the underlining group is 
abelian, the weight M aia2a3 (z) is automatically cyclic and gauge invariant. Further- 
more, M aia2a3 (z) does not depend on the orientation of Fk- 

The choice of z in (fT6|) will determine the model we are describing. For example, 
consider the curve (a , a 1 ) = (| e j \e~^) parametrized by /3. The corresponding 
weight 

M aia2ai {\eP ,\e-?) = e «-D— (17) 

corresponds to the spin-gauge model. We can also choose the curve (a , a 1 ) = 
(2' \ e ~ 2 ^) an( ^ that wm S^ ve 

M aia2a3 (^e-^) = e-^-(-D- + — »], ( 18) 

describing yet another model. 

The second algebraic information is a co-product A : A — > A <%) A defined by 

A(0 C ) = A c ab a ® 4. (19) 

The tensor A c afe defines also a product on the dual space A*. Using the dual basis 
{ip a } we define 

ifj a ■ V b = A c a V- (20) 
In analogy with (fT5|) and (|16p we define the co-trace T* £ ^4 as 

T* = A b fc > a . (21) 

and the tensor A bl "' bn as 

A 6i-6„ _ ^01 .^...^y*). (22) 
It turns out that the co-product that we will need is very simple. We will set 

A c ab = 5 a J b c . (23) 

Therefore 

£bi-b n = [ 1 if a11 indices are equal 



otherwise 

Note that the orientations of hinges do not affect the corresponding statistical weight. 

The tensors M a b c and A a bc given in (|13p and (f23|) together with the weights M a f, c 
and A ai "' an defined by (|16p and (|22p completely specify our model. It is a simple 
exercise to show that the partition function ([9]) reduces to 

Z(a°,a 1 ) = J2 II M ^a 3 (f), (25) 
fo}/eT(M) 

where the sum {07} is over the configurations on the links / 6 T(.M) and M a b c (f) is 
the weight of a configuration (01,02,03) on the face / £ T(.M). Notice that 

^(Aa^Aa 1 ) = A^Z(a°, a 1 ). (26) 
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As M aia2a3 depends on z = a°(j)Q + a l 4>i G A, the model depends on two param- 
eters (a , a 1 ). We will denote the parameter space by V . The parameters (a , a 1 ) 
are related to the parameters (mo, mi) of ([7]) as 

too = a + a 1 

mi = a ° -a 1 . (27) 

Let us recall that the algebra defined in (|12p is a group algebra and as such 
it is also a Hopf algebra with co- product coming from A a bc . The maps antipode 
S : A — > A, unity e : C — > A and co-unity e : A — > C can be described in terms of 
tensors as S((j) a ) = S b 4>b, e(l) = e a (/) a and e((j) a ) = e a . For the case of Z2 we have 

S b = <t e a = 5° a and e a = 1. (28) 

Notice that these tensors are essentially trivial and will not show up explicitly in the 
calculations. For a non abelian case, for example, S b is related to the orientation 
but that will play no role in the Z2 case. 



3 Quasi-Topological Limits 

We would like to explore the model (|25p and look for points of the parameter space 
(a , a 1 ) where the model has a topological or quasi-topological behaviour. The 
simplest case is the one considered in [HE]. It corresponds the point (a , a 1 ) = (1, 0) 
in the parameter space T or, equivalently, to the choice in (fT6|) of z equals to the 
identity (fio of the algebra. For this particular point of the parameter space we can 
bring in the Hopf structure of A, follow [JJ [5] and conclude that 

Z (T(M), a = 1, a 1 = 0) = 2 Nf+Nl ~ Nt Z top (M) , (29) 

where Nt, Nf and Nl are the number of tetrahedra, faces and links of T(A4). 
However, (a , a 1 ) = (1,0) is not the only quasi-topological point. In this section we 
will show that there is a sub-set To C T with dimension one such that the partition 
function is quasi-topological. In figure [1] we have the parameter space V where the 
set To consists of four straight lines: the diagonals and the axis. We also have 
included the models with weights W^- 1 ' and defined in ^) and j5]). They are 

curves parametrized by a single parameter /3. 

Let us consider (a , a 1 ) = (A,0). It corresponds to z = Xcpo G A. The new 
weight is simple M a i, c {z) = AM a b c (0o). Therefore the partition function is the same 
as for (a , a 1 ) = (1,0) multiplied by a factor. In other words 

Z(T(M), A,0) = \NF 2 NF+N L -N TZ top _ ( 30 ) 

The tensor components M a \, c {z) and A bl "' bn depend on the choice of a basis 
of the algebra A. Therefore, different choices of basis will lead to different weights 
and therefore different models. However, the partition function Z(T(M),aP , a 1 ) has 
been written as a scalar and therefore is invariant under a change of basis. This large 
invariance of Z(T(A4), a , a 1 ) has been interpreted as dualities between different 
models. This fact has been explored by us in |11] to show that the classical Kramers 
and Wannier duality relations are special cases of these more general dualities. In 
what follows we will show that there is a duality relation between the model at 
(a , a 1 ) = (A,0) and (a , a 1 ) = (0,A). We will show that 

Z(T(M),0,X) = Z(T(M),X,0) , (31) 



8 




Figure 4: The set T consists of the two axes plus the dashed diagonal lines. The models 
and are curves on T. 



which allow us to compute Z(T(A4), 0, A). 

Let us start by recalling that at the topological point (a , a 1 ) = (1,0) the weights 
read 

M ahc {<Pv) = 25{a + b + c,0) , (32) 
A bl "' bn = 5(b 1 ,b n )5(b 2 ,b n )---5(b n . u b n ). (33) 



Consider the matrix 



L-[ I I ) . (31) 



and a new basis {4>' , 4>i} defined as <j>' a = E^cpb- This transformation simply changes 
an index a by a, where 0=1 and 1=0. In the new basis we have 

M' abc (l, 0) = M s5g (l, 0) = 25(a + b + c,0) = 25(a + b + c, 1) , (35) 

= 6(b 1 ,b n )5(b2,b n )---8(b n -i,b n ) . (36) 

Note that 

£\jbi---b n _ ^bi—b„ (37) 

On the other hand, the tensor components M a b c (0, 1) on the original basis is M a b c (0, 1) 
S(a + b + c, 1). Therefore 

M^ C (1,0) =M afec (0,l). (38) 

Equations ([371) and (EHJ) imply that Z(T(M),0,X) = Z(T(M),X,0) as announced. 

The partition function corresponding to the diagonal line a = a 1 in figured] can 
be obtained by choosing Zi n t = (fro + <f>i in (|16p . One can see that Zi n t is such that 
0g 'Zint = Zi n t- In a Hopf algebra, such element is called a co-integral [12] . Therefore 

M abc (l,l) = (T,z in t) = 2. (39) 

The fact that the weight M abc {l, 1) is independent of the configurations turns the 
computation of Z(T(A4), A, A) completely trivial. One can immediately see that 

Z(T(M),\,X) = 2 Nf+Nl X Nf . (40) 
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Note that Z(T(A4), A, A) depends only on the size of the lattice T(M). In contrast 
with (j30|) . the topological invariant is trivial since the partition function does not to 
depend on the topology of the underlining manifold M. 

To compute the partition function for a 1 = —a it is enough to consider the 
point (a , a 1 ) = (1,-1). That is the same as setting z = <pQ — §\ on (fT6|) . The 
weight M a b c (l,— 1) reads 

M afec (l,-l)=2(-ir+ 6+c . (41) 

Instead of computing Z(T(A4), 1, —1) directly, we will make use of the duality related 
to change of basis. Let us choose another basis by applying the transformation 
matrix E h a = (— l) a S\. In the new basis the weights are 

M' abc = 2, 

{ ^\J\aia2---a n _ / _-j^ai+a 2 H ha n jai fid2 . . . fiQ>n-l (42) 

The resulting model has weights associated to the links only and the partition func- 
tion can be easily computed. After plugging (j42j) in ([9]) and taking into account the 
scaling factor we have 

Z(T(M), A, -A) = A^ + C- 1 )^ ) > ( 43 ) 

i 

where the product runs over the links of the lattice and 1(1) denotes the number of 
faces that share the link I. Note that 1(1) is not of topological nature. Furthermore, 
the partition function vanishes whenever there is a link that is shared by a odd 
number of faces. This is an indication that (a , a 1 ) = (A, —A) is a very peculiar 
model. 

It is clear from the computation of the partition function that the lines that make 
up Tq are not all equivalent. Actually they are all different from each other. It is 
true that Z(T(A4), A, 0) is the same as Z(T(A4), 0, A). However, the expectation 
value of Wilson loops are not the same for these two models. As we will show in the 
next section, only (A, 0) is in fact a topological theory. 

Before we conclude this section, we would like to point out the relation between 
To and familiar models such as the ones defined by (|5J (spin-gauge model) and © 
(gauge theories). As we have discussed before, the spin-gauge model given in (fT7|) 
corresponds to the curve (a , a 1 ) = (l/^e' 3 , l/^e - ' 3 ). It is clear from figure [4] that 
this curve approaches Tq as (3 goes to +oo and — oo. Another point of contact with 
Tq is when /3 — > 0. As for the model ©, only the limits (3 — > and (3 —> ±oo are 
part of Tq. 

4 Wilson Loops 

The two parameter gauge model (|10p of last section have numerical quantities 
(Wit (7)) that are the natural generalization of the expectation value of Wilson loops 
for a closed curve 7 and irreducible representation R. The definition of (Wr(7)) re- 
duces to the familiar expression when restricted to the usual gauge models. In this 
section we will define and compute these observables for To. 

For simplicity, we will start by considering the loop 7 to be unknotted. Knotted 
loops will be considered in the last part of this section. 

Let 7 be a loop made of a set of links uj(j) = {u>i, ...,U) p }. For such a loop we 
introduce the tensor W aia2 ... ap with p indices given by 

W*.. ap = (T,<f> ai ...<t> ap -z R ), (44) 
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where is the configuration at link Ik and zr is the unique element in the center 
of A such that 

(T, 4> g ■ zr) = Tift(g), (45) 

where R denotes an irreducible representation of the group. 

The group Z2 has only two irreducible representations labelled R = and R = 1, 
such that 

Tr (^a) = 1 , 
Th(^) = (-l)V 

We only need to consider the non-trivial representation. In other words, we will set 
zr = 7; (4>o ~ <t>i) since that will give us (T,(f) a ■ zr) = ( — l) a . For now on we will 
omit the index R indicating the representation simply write (|44p as 

W ai ... ap = (_l) a *+~+ a P. (46) 

We would like to construct (W(j)) as a scalar in the same way it has been done 
for the partition function (fTUj) . As before, we make use of the contra- variant tensors 
A bl " bn associated to the links. Consider a link I shared by 1(1) faces. If I does 
not belong to the loop 7, the corresponding tensor is the same as for the partition 
function and will be written as /\ bl "' b iW . If, however, the link in question is one of 
the links of 7 (I = ujj € ^(7)), the corresponding tensor will be A fcl '" fcl W aj '. Note 
that the new tensor has an extra contra- variant index a,-. The expectation value of 
the Wilson loop is defined to be 

(^(7)>=§ n m ^(/) n ^"Mo n ^■■■ b ^ Wav .. ap , m 

where we are using the simplified notation of (|10p . The contraction of indices are 
as follows: each covariant index of W ai ... ap is to be contracted with the extra index 
dj in A bl '" bl( - l '> a: > as explicitly indicated. The indexes from the product ri/e{F} M a bc 
follows the same rule as for the partition function when contracting with the bk 
indexes in A bl '" b i(i) and A bl '" b i{i) a i . An interpretation of ([4T|) in terms of gluing of 
hinges (H m ) and faces (Fk) that is analogue to the partition function can also be 
given. To each link I we associate a hinge Hi. If the link is not part of the loop 7, the 
corresponding hinge has exactly 1(1) flaps that will connect to 1(1) faces in the usual 
way (see figure [3|). For links ujj £ ^(7) the corresponding hinge H Uj has (I(uJj) + 1) 
flaps. After gluing the I(u>j) faces to H u ., we are left with extra flaps, one for each 
link of a; (7) as illustrated in figure [5j The contraction of indices a,j in (|47p can be 
seen as the attachment of a polyhedral face W 7 with p edges. An example of such 
attachment is shown in figure [6j This special face W 7 has weight W ai ... a and is not 
to be thought as part of the lattice. In general, it will not be possible to embed 
in 3D space. 

The expectation value of the Wilson loop will be a function (W('y))(a , a 1 ). As 
a consequence of (|26p one can see that 

<W(7))(Aa°, Aa 1 ) = {W(i))(a°, a 1 ). (48) 

Therefore we only need to compute (W( , y))(a , a 1 ) for one point on each straight 
line of Tq . 

We now proceed with the computation of (W(7)) when the parameters (a , a 1 ) 
belongs to Tq. That can be divided in three cases as follows. 
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(a) (b) 

Figure 5: A small region of the lattice with a Wilson loop 7 is illustrated in (a). Figure 
(b) shows the corresponding decomposition. For a link of 7, the corresponding hinge has 
an extra flap as indicated. 





Figure 6: Wilson loops can be thought as an extra face which is not part of the lattice. 



Case 1 (z = (fto + <pi): 

The first case correspond to (a , a 1 ) = (1, 1). As we have seen, M a b c = 2. All the 
sums on b indices in ([47p are straightforward due to the Kronecker deltas in A r " J (0 
and A 6l '" 6l «°J. We are then left with 

1=1 Oj 

This result does not depend on the loop 7. 

Case 2 (z = 4>o — <j)\): 

For this particular case we have Wai-op = (-l) ai +-+V Notice that this is the 
same function as the weights M a t, c = 2( — l) a+b+c for the faces. Therefore the numer- 
ator in (T47J) is the same thing as a partition function with an extra face determined 
by the loop 7 (see figure [6]) . Using (|43p we can write 

(wm = m — \ 1)m n [i + n [i + • 
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This region of the parameter space is quite peculiar. Notice that the denominator 
of (W(7)) vanish if /(/) is odd for some I and (W("f)) is not well defined. When 1(1) 
is even for all links, (W(j)) is well defined but it is equal to zero due to the factors 
[1 + (_!)/«+!]. 

Case 3 (z = <ft g , g = 0, 1): 

The points with coordinates (a , a 1 ) given by (1,0) and (0,1) in Tq correspond 
to z = (j) g with g = and g = 1 respectively. We know from previews sections that 
the partition function is the same for these two cases. However, this is not true for 
(W(j)). 

We will show that (W(j)) is quasi-topological in the sense that it does not 
depend strongly on the geometry of 7. The idea is to investigate the behaviour 
of (W(y)) under small deformations. In a triangulation it is natural to define a 
small deformation of a loop 7 as follows. Let 01(7) be the set of links of 7. A 
small deformation of 7 is a local move that replaces one link u)\ € ^(7) by a pair 
of links W2>W3 such that 0J\,0J2 and U3 belong to the same face. There is also the 
reverse move where a pair of links uj\,U2 is replaced by W3 provided they belong to 
the same face. A small deformation is shown in figure [7J With this definition of a 
small deformation we now establish how (W(^)) changes under small deformations 
for z = <j) g ,g = 0, 1. 




(a) (b) 

Figure 7: The loop 7' in (b) is a small deformation of the loop 7 in (a). 



Consider figure [8] (a) where we have single out a small part of the partition 
(Fk,Hj). Only the elements connected to the link u\ are relevant. In figure [S] (b) 
we have performed a small deformation of 7 by replacing L02, W3 by oj\. The weights 
associated to figure [8] (c) and figure [8] (d) are the tensors A and B given by 

,a2 - a v b2---b q C2---c r _ ,^ » ai—a p u>i \ bv-bq A ci—CrXir 

^ai-as — i "ai&ici^ A yv u>\ai— a s 

= M aiblCl A"?A b ^A c ^W aiai ... as (51) 

ai,&i,ci 

and 

Tja2—a p b2—b q C2—Cr _ A ai A 61 ■• '6 4 W2 A Ci -Cr4)3 TI/ 

= Y, M ai6lCl A a i- a fA 6l - 6 «A Cl -^W 6lCiai ... Qs , (52) 

ai,bi,ci 
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(d) 

Figure 8: In (a)-(b) we represent a small deformation of a loop 7 in a lattice and in 
(c)-(d) we represent the same deformation in a lattice decomposition. 



where we have performed the sums on uj±,U2 and W3. We also have explicitly written 
the sums on ai,&i and c\. Note that W aiai ... as = W ai ... as , W biciai ... as = 

(— l) 6l+Cl W Q1 ... Qs and M ai 6 lCl = 25{a\ + b\ + ci, That is enough to show that 

M aiblCl W bl = (-l) 9 M aiblCl W aiai ... as , (53) 

therefore 

Aa2---a p b2---b q C2---c r _ , -,\g Tj a 2---a v b2---b q C2---c r fKA\ 

S± ai ... as — { — !) r>ai—a a ■ \ oq: ) 

We see that {W(j)) is invariant under small deformations for g = that corre- 
spond to (a , a 1 ) = (1,0). As for g = 1 or (a , a 1 ) = (0,1), the number (W^)} 
flips sign each time we perform a small deformation on 7. 

The simplest loop 70 is made of a single triangular face with links (^(7) = 
{wi,W2,w 3 }. If we recall that M ahc = 25(a + b + c,g) and W mW2UJ3 = (_i)^i+^2+^ 3 
, it becomes a straightforward computation to show that (II / (7o)) = ( — l) 9 - We can 
now deform the smallest loop 70 by adding (N — 1) triangles and arriving at a planar 
loop 7at. For such a loop we get 

(W( 7N )) = (-iy N ; 5 = 0,1. (55) 
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Equation (|55j) for g = 1 shows (W(7jv)) as a function of the number of triangles 
swept in the process of stretching 70 into 7jy. This function is a very simple "area 
law". It depends only on the parity of N. Notice that each time we add a triangle 
to 70, the number of links of the loop also changes by one unity. Therefore we could 
have written 

(W( 1N )) = (-1)^ , (56) 

where Nl is the number of links of the loop 7. We could interpret this formula as a 
"perimeter law". The fact that there is no distinction between area or perimeter law 
is peculiar to the gauge group Z2 and the fact that we are using triangular lattices. 
For square lattices the variation in the number of links is even and (|56|1 does not 
hold, but (|55p is still true. 

So far we have considered only planar loops. For knotted loops, the expectation 
value (W(7jv)) f° r 9 = may depends on the class of isotopy of the loop 7. Equation 
(|55|) is valid for loops 7 that can be deformed to the trivial knot. If 7 is a non trivial 
knot, equation ([55]) may not hold. We will show that (W(7jv)) actually does not 
depend on the class of isotopy of 7 and therefore (|55p is still correct. 

Let 7 be a knotted loop and its knot diagram as illustrated in figure [9] (a) . We 
can produce a new knot diagram by flipping under-crossing into over-crossings and 
vice versa. This flips are local moves that only affect the knot in a small region. It 
is a well known result from knot theory that any 7 can be made into the trivial knot 
if we perform a number of flips, as we can see in figure [9] (b). We will show that 
(W(7jv)) is invariant by flips and therefore does not depend on the isotopy class of 
7- 




(b) 

Figure 9: In (a) we have a knotted loop that can be made trivial by two flip moves, as 
illustrated in (b). 

Let us consider a 3-dimensional ball B around a crossing in a knot K. That 
will give us curves 7^ and 7^ connecting the points (a, b) and (c, d) at the surface 
dB = S 2 of B as in figure [TOl (a) , After a flip move, we have a new knot K and new 
curves are 7^ and %d as illustrated in [10] (b). Before analysing the general case, 
we will look at a simple example where B is a cube inside the triangulation and the 
curves, before and after the flip move, are the ones given in figure [TT1 (a) and (b). 
In this figures we are using the interpretation of the expectation value of the Wilson 
loop as an extra flap as was explained in section [4[ 
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d b d b 




a c a c 



(a) (b) 

Figure 10: The situation shown in (a) represents a flip move of one shown in (b) and 
vice versa. 




One can see that the the sequence of extra flaps can be deformed into the se- 
quences given by figure [121 (a) and (b). After we perform these deformations, the 
computation of (W(K)) and of (W(K)) differ only at a single link. The difference 
is that the flaps coming from the two curves are swaped. In one case, the curves will 
contribute to (|47p with a factor 

A vv u 1 a 1 ---a s VV UJ2/3l ... l3r , [01) 

and in the other case the factor is 

^ VVui^-asVV^^...^- (7°J 

The tensor /\ a i ,, ' a ""'2Wi j s invariant by permutation of the indices uj\ and 0J2 and 
these two factors are the same. Therefore, when comparing (W(K)) and (W(K)) 
we can use a sequence of small deformations to arrive at the configuration on figure 
[TUl (a) and figure [TUl (b). The flip itself will not give any contribution and we can 
use the result for planar loops given in (|55j) . 

As for the generic case we can proceed as follows. Choose curves 71,72 and 73 
connecting the pairs (a, d), (d, b) and (c, b) as in figure [T3l Since the 3-ball is simple 
connected, 7^ is isotopic to 71 072 and 7^ can be deformed to 73 . In a similar 
way, we have ~{ a b is isotopic to 71 o 72 and 7^ can be deformed to 73 o 7^ . In a 
similar fashion as for the particular case of figure [10] (a) , there will be two sequence 
of flaps along 72. One comes from the the deformation of 7^ and the other comes 
from the deformation That has to be compared with a similar sequence of flaps 
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(a) 



(b) 




Figure 12: The pictures (a) and (b) differs from each other only by the link (d, b), as we 
can see in (c) and (d). 




Figure 13: The paths 7 1; 72 and 73 are some paths that connect the vertices (a, d), (d, b) 
and (b, c). 



coming from the deformation of 'jab and 7 cc z- As in the case of figure [T2l (c) and (d), 
these sequences of pairs of flaps along 72 can only differ by a permutation. Since 
A" 1 '""" is symmetric by permutation of indices we can conclude that flips do not 
give any contribution and we can use the result for the planar loops given in (j55|) 
for any knot K. 
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5 Final Remarks 



Topological field theories are among the simplest lattice models we can have. From 
the physics point of view they are peculiar models. Partition function and corre- 
lations can be computed but the dynamics is too simple. Rather than considering 
TQFTs in isolation, we have looked at the problem from a broad perspective and 
investigated a two parameter family of models where TQFTs can arise at certain 
points of the parameter space. We have considered gauge theories with Z2 sym- 
metry since it is the simplest gauge group but can still accommodate non trivial 
models, such as the 3D spin-gauge model. These more familiar models appear as 
one parameter curves in the two dimensional parameter space V. 

We have found several limits that we can loosely call topological or quasi-topological 
comprising a subset Tq of T. On Tq both partition function and expectation value 
of Wilson loops were computed. The partition function points on Tq are topological 
numbers up to an overall scale factor. One could think that Tq contain only topo- 
logical models but the expectation value of the Wilson reveals something else. First 
of all, (W i?(7)) does not depend on the isotopy class of the curve 7. Furthermore, 
for a subset of To, {Wr(^)) depends on the size of 7 and follows a discrete version 
of an area law. 

In the parametrization (a , a 1 ) of V used in the paper, the subset Tq is made 
of 4 straight lines passing through (0,0). By looking at the gauge Ising model, we 
can see that it approaches three of these lines for (3 — > ±00 and f3 — > 0. There is an 
extra line given by a 1 = —a that, as far as we know, does not relate directly to 
any physical model. 

The existence of a set Tq in the parameter space where the model behaves in a 
topological way can be seen as an Euclidean version of topological order. It seems 
that, rather than a special case, the same phenomena will happen for gauge theories 
with any compact gauge group G. For Zjy and non abelian groups, the analysis is 
much more involved and it will be reported in a separated paper. 
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